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Abstract. We study reductions of the Hamiltonian flows restricted to their 
invariant submanifolds. As examples, we consider partial Lagrange-Routh 
reductions of the natural mechanical systems such as geodesic flows on compact 
Lie groups and n-dimensional variants of the classical Hess-Appel'rot case of 
a heavy rigid body motion about a fixed point. 
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1. Introduction 

In this paper we study reductions of the Hamiltonian flows restricted to their 
invariant submanifolds. Apparently, the lowering of order in Hamiltonian systems 
having invariant relations was firstly studied by Levi-Civita (e.g., see |17j . ch. X). 



1.1. Hess-Appel'rot System. The classical example of the system having an 
invariant relation is a celebrated Hess-Appel'rot case of a heavy rigid body motion 

1 



2 



BOZIDAR JOVANOVIC 



[T11[T]. Recall that the motion of a heavy rigid body around a fixed point, in the 
moving frame, is represented by the Euler-Poisson equations 

(1) — rfi = m X + (39Jt7 x r, — 7 — ^ x Co, uj — Am, 
at at 

where uj is the angular velocity, m the angular momentum, I — A^^ the inertial 
tensor, QJt mass and r the vector of the mass center of a rigid body; 7 is the direction 
of the homogeneous gravitational field and <S is the gravitational constant. 

The equations ^ always have three integrals, the energy, geometric integral and 
the projection of angular momentum: 

(2) - i(77i, Lu) + m&ir, 7), ^2 = (7, 7) = 1, -^3 = {m, 7). 

For the integrability we need a forth integral. There are three famous integrable 
cases: Euler, Lagrange and Kowalevskaya plfM]. 

Apart of these cases, there are various particular solutions (e.g., see [11]). The 
celebrated is partially integrable Hess-Appel'rot case [mH]. Under the conditions: 



(3) ^2 = 0, riy/as - 02 ± r3y/a2 - ai = 0, A = diag(ai, 02, 03), 
a^ > a2 > ai > 0, system ([T]) has an invariant relation given by 

(4) = (m, r) = rimi + r^ms = 0. 

The system is integrable up to one quadrature: the compact connected compo- 
nents of the regular invariant sets J-'i = c\,T2 = 1,^3 = 23,^-4 = are tori, but 
not with quasi-periodic dynamics. The classical and algebro-geometric integration 
can be found in jT4] and [9], respectively. 

There is a nice geometrical interpretation of the conditions ([3]): the intersection 
of the plane orthogonal to r with the ellipsoid (Am, m) = const is a circle (e.g., see 
[8]). If instead of the moving base given by the main axes of the inertia, we take the 
moving base /i , /2 , fs such that the mass center of a rigid body f is proportional 
to fs, 



(5) r^pfs, p^yrl + rl, 
then the inverse of the inertial operator reads 

f 0.2 ai3 

(6) A= a2 023 

\ ai3 123 133 

and the invariant relation ([4]) is simply given by 

(7) 7713 = 0. 

Historical overview, with an application of Levi-Civita ideas to the Hess-Appel'rot 
system and other classical rigid body problems can be found in Borisov and Mamaev 
0. 
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1.2. Reduction of Symmetries. Let G be a Lie group with a free proper Hamil- 
tonian action on a symplcctic manifold (Af, cj). Let 

(8) $ : Af 0* 

be the corresponding equivariant momentum map. Assume that i] is a regular value 
of so that 

(9) M, = 

and Mq^ = $"^(0^) are smooth manifolds. Here Oj^ = G/Gn is the coadjoint orbit 
of The manifolds and Mq^ are Gj^-invariant and G-invariant, respectively. 
There is a unique symplectic structure w,, on Nrj = MrJGr^ satisfying 

(10) uj\m„ = d7r*a;,, 

where tt^ : N^i is the natural projection f20]. 

According to Noether's theorem, if is a G-invariant function then the momen- 
tum mapping $ is an integral of the Hamiltonian system 

(11) x^Xh. 

In addition, its restriction to the invariant submanifold Mri projects to the Hamil- 
tonian system y — Xh on the reduced space Nn with H defined by [50] 

(12) h\M,^H oTi^. 

An alternative description of the reduced space is as follows. Let {•,•} be the 
canonical Poisson bracket on (M, w). Then the manifold M/G carries the induced 
Poisson structure and Mq^/G is the symplectic leaf in M/G. The mapping 

(13) * : ^ MoJG 

which assigns to the G,,-orbit of x G the G-orbit through x in Mo^ establish 
the symplectomorphism between A'^ and Mq^/G (e.g., see [21])- 

1.3. Outline and Results of the Paper. The aim of this paper is to study 
Hamiltonian systems that naturally generalize geometrical properties of the Hess- 
Appel'rot system (Definition [31 Section 4). This is the reason we are interested 
in the following: suppose that h is not a G-invariant function, but ([Hj) is still an 
invariant manifold of the Hamiltonian system (jlip . As a modification of the regular 
Marsden-Weinstein reduction, we study (partial) reduction of the Hamiltonian sys- 
tem (fll]) from Mri to Nr, and relationship between the integrability of the reduced 
and nonreduced system (Theorems [TJ [21 Section 2) . 

In Section 3 we define partial Lagrange-Routh reductions (Theorem [3]). The 
construction of the geodesic flows on compact Lie groups, having invariant manifolds 
foliated by invariant tori is presented. However, the flows over invariant tori are 
not quasi-periodic. 

In Section 4 we consider an n-dimensional heavy rigid body motion about a fixed 
point. Recall that there are two natural multidimensional variants of a heavy rigid 
body which are described by the Euler-Poisson equations on the dual spaces of 
semi-direct products so{n) x so{n) and so(n) x R". The integrable generalizations 
of the Lagrange top with respect to the first and the second n-dimensional variant 
are given in [25|, 1 10] and [3] , respectively. The n-dimensional Kowalevskaya top in 
the presence of q gravitational fields {q < n) represented by Euler-Poisson equations 
on the dual space of the semi-direct product so(n) x (R")' is given in [26] . 
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Recently, an interesting construction of the Hess-AppeProt system on so{n) x 
so{n) is studied by Dragovic and Gajic Besides, by generalizing analytical 

and algebraic properties of the classical system ([T|), Q, they introduced a class of 
systems with invariant relations, so called Hess-AppeVrot type systems with remark- 
able property: there exist a pair of compatible Poisson structures, such that the 
system is Hamiltonian with respect to the first structure and invariant relations are 
Casimir functions with respect to the second structure (for more details see 

We consider the another generalization of a heavy rigid body and define the 
n-dimensional Hess-Appel'rot system within the framework of partial reductions 
(Lemma |6]). It is known that the motion of a rigid body mass center in the classical 
Hess-Appel'rot system is described by the spherical pendulum equations [29l [8]. 
We shall prove the same statement for the n-dimensional variant of the system 
(Theorem 111)- On the other side, the L-A pair of the system written in the closed 
form of the Euler-Poisson equations on the dual space of the semi-direct product 
so{n) X M" is given. In particular, we get the L-A pair representation for the Euler- 
Poisson equations of the classical Hess-Appel'rot system ([I}, (O, ([6]) different from 
those given in [9] (Theorem [9|) . 

Finally, we have made a relationship between our construction and the Dragovic- 
Gajic so(4) x so(4) system (Theorem [TUl Section 5). 

2. Partial Reductions and Integrability 

Let G be a compact connected Lie group with a free Hamiltonian action on 
a symplectic manifold (M, w) with the momentum map ([8]). Assume that 77 is a 
regular value of <&. 

Theorem 1. (i) Suppose that the restriction ofh to Mq^ is a G -invariant function. 
Then Mr/ is an invariant manifold of the Hamiltonian system ill}) and X^Im^ 
projects to the Hamiltonian vector field Xh 

(14) d-n:y^{Xh)\x = XH\y=-,r^{x)i 

where H is the induced function on defined by il2\) . 

(ii) The inverse statement also holds: if (0) is an invariant submanifold of the 
Hamiltonian system 177]) then the restriction ofh to M^i is a Gjj-invariant function 
and Xh\Mr, projects to the Hamiltonian vector field Xh on N^j, where H is defined 

by m- 

In both cases, the Hamiltonian vector field Xh is not assumed to be G-invariant 
on M. Moreover Xh\M may not be G^-invariant as well. It is invariant modulo 
the kernel of duri, which is sufficient the tools of symplectic reduction are still 
applicable. 

Definition 1. We shall refer to the passing from x = Xh\Mn to 

(15) y = XH 
as a partial reduction. 

If h is G-invariant, the complete integrability of the reduced and original system 
are closely related (see [16j[28]). In our case, we have the following corollary. 

Theorem 2. (i) Suppose that the partially reduced system is completely in- 
tegrable, with a complete set of commuting integrals Fi,...,Fm, m = ^dimiV^. 
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Let fi = Fi o TT^, . . . , /„ = Fn o T^ri- Then Mri is almost everywhere foliated by 
(m + dim Gri)- dimensional invariant isotropic manifolds 

(16) Mc = {/i = ci, . . . , /,„ = c„i} 
of the system ill]) . 

(ii) In addition, if rj is a regular element of g* and /i, . . . , can be extended to 
the commuting G -invariant functions in some neighborhood ofM^i, then the compact 
connected components of invariant manifolds Hid]) are tori. However, in general, 
the flow over the tori is not quasi-periodic. 

Remark 1 . If the partially reduced system (jlSp is integrable then we need (dim G^)- 
additional quadratures for the solving of i = Xhlm^ {the reconstruction equations). 
Also, if Gjj = G, i.e., O,, = {77}, then m+dim = i dimM. Whence, the invariant 
manifolds (|16p are Lagrangian. 

The partial reduction can be seen as a special case of the symplectic reductions 
studied in [Ulin] (see also [TB], Ch. III). There, the symplectic reduction of a sym- 
plectic manifold (M, w) is any surjective submersion p : iV — > P of a submanifold 
N <Z M onto another symplectic manifold (P, f2), which satisfies p*Vt = aj|jv- Also, 
instead of a compact group action one can consider a proper group action. However 
we work in the framework of a regular Marsden-Weinstein reduction and a compact 
group action, which allows us to easily describe the partial reductions of natural 
mechanical systems considered in this paper. 

Proof of Theorem\^ (i) Let ^1, . . . , be the base of g, (77, ^q.) = r^a and 

(17) = ($,^a) : A/^M, a = 

Then the level set of the momentum mapping ([9]) is given by the equations 

= ?7q, a = 1, . . . ,n. 

The action of G is generated by Hamiltonian vector fields X^^. Since h\Mo is 
G- invariant, for x S we have 

(18) [dh,X^J^{h,(t)^}^-{d(t)c.,Xh)^Q, a = l,...,n. 

Thus Mr, is an invariant submanifold. 

Let h* be an arbitrary G-invariant function that coincides with h on Mq . Then 
Xh- is a G,,-invariant vector field on My, which project to Xh [20] : 

(19) dny,{Xh')\x = XH\y=T,^[x)- 

Let S = h — h* . From the condition (5|j\/^ = we can express S, in a an open 
neighborhood of Mj,, as 

n 

a=l 

Now, let / be a G-invariant function on M. Since {/,0a} 
theorem), we get 



= (Noether's 
= 0. 
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Thus the Poisson bracket {(5, /}|a/^ = ^{df, Xs)\Mr, vanish for an arbitrary G- 
invariant function /. In other words 

(20) XsUeT4G-x)nT^Mr,. 

Combining (fT9|) . Xh = Xw + Xg and (f20|) with the weU known identity (e.g, see 

M) 

Tx{G ■ x) n T.j:Mri = Tx{G,j ■ x) = kerdTT^Ia; , 
we prove the relation (I14p . 

(ii) If G is a connected compact group, the coadjoint isotropy group is 
connected. Thus, the function h is G^-invariant if and only if it is invariant with 
respect to the infinitesimal action of G^. 

Suppose that Af^ is an invariant submanifold of (fTT|) . Then holds for x e M^. 
Since the action of G is generated by Hamiltonian vector fields X^,^ , from ^TE\\ we 
get that h is invariant with repsect to the infinitesimal action of G,,. Therefore we 
have well defined reduced Hamiltonian function H on N^. 

By using the diffeomorphism (fT3|) and the fact that Mo^ is a closed submanifold 
of M, we can find a G-invariant function h* on M which coincides with h on M^. 
Now, the relation (|14p follows from the proof of item (i) . □ 

Proof of Theorem [H (i) Consider a regular invariant Lagrangian submanifold 

(21) ^f^ = {Fl=Cl,...,Frn= C„} C Nr,, 

of the partially reduced system (fT5)) . From the relations (jlOp and we get that 
A4c = ""ly^lM;) is an invariant isotropic manifold of the system (fTTI) . 

(ii) Let 77 e 0* be a regular element (G,, « T'' is a maximal torus of G, 
r = rank G). If the connected component Af° of (PT|) is compact then, by Liouville's 
theorem, it is diffeomorphic to a m-dimensional torus T"* with quasi-periodic flow 
of (fTS]) . Thus the compact connected component M° = n^^{N'°) of (fTH]) is a torus 
bundle over T™: 

T'' — > 

(22) i TTr, 

Let Ji , . . . , Jy, be the basic Adg-invariant polynomials on q* and Pi , . . . , Pn-r be 
linear functions on g* such that Pfc, /q are independent at rj. Also, let = /„ o <I> 
and pfe = Pfc o $ be the pull-backs of and Pfe by the momentum mapping, 
a = l,...,r, k = 1, . . . ,n — r. 

Suppose that /i , . . . , /,„ can be extended to commuting G-invariant functions 
in some G-invariant neighborhood V of M°. Then, within V, Ai° is given by the 
equations 

fi = ci, ...,/,„ = c„i, ii = 11(77), . . . , V = Irii]), pi = Pi (77), . . . ,P„_4(?7) ■ 

From the Noether theorem the functions ia, Pk commute with all G-invariant 
functions on M and, since ia are G-invariant, the following commuting relations 
hold on V: 

(23) {fajb} = {fa,ia} = {«a,«/3} = {fa,Pk} = {ia,Pk} = 0, 

a, 6 = 1, . . . , TO, a, /? = 1, . . . , r, k — l,...,n — r. 

Now, as in the case of non-commutative integrability of Hamiltonian systems 
[131 [53], M° is a torus with tangent space spanned hy Xf^, Xi^ . Namely, from the 
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vanishing of the Poisson brackets ([23|) we get that the vector fields Xf^, Xi^ are 
tangential to M°. Since they are independent, from the dimensional reasons, they 
span the tangent spaces T^M", x € A4°. Taking into account and the relations 
Lo{Xf,Xg) = [Xf,Xg] = Xy^g^, f,g& C^(M), we (re)obtain that 

is isotropic. Furthermore the vector fields Xi^ commute between themselves. 
Since is a compact manifold admitting rn+r = dim independent commuting 
vector fields it is a (m + r)-dimensional torus, i.e., the bundle ((22|l is trivial. 

In general, the flow of i = Xh over the torus Al° is not quasi-periodic: the vector 
field Xh do not commute with vector fields , . . . , Xf^ , Xi-^ , ■ ■ ■ , Xi^ (although 
Poisson brackets {h, fa}, {h, ia} vanish on □ 

Remark 2. Let x e A4° and z = tt{x), where tt : M ^ M/G is the canonical 
projection. By the use of ([13]) and the existence of local canonical coordinates on 
the Poisson manifold M/G within some neighborhood U oi z (e.g, see [H]), the 
functions fi can be always extended to G-invariant commuting functions in the 
G-invariant neighborhood V = ti^^{U). Thus, if M° is "small enough", it is a 
torus. 

Remark 3. It is clear that complete commutative integrability of the reduced system 
in Theorem [5J can be replaced by the condition of non-commutative integrabilty. 

Remark 4. The Hamiltonian flow preserves the canonical measure — 
(the Liouville theorem). If A/^ is an invariant manifold of the equations (jlip . then 
the functions (jl7[) are particular integrals. Therefore the time derivative of 0q is of 
the form 

(t)a{x) = h} 1paf3ix){(l)a - Ca), 

P 

where ipaf3 are smooth functions. Let trip be the trace of the matrix ipaf3- After 
straightforward calculations, one can prove that the flow pT|) restricted to the 
invariant manifold Af^ preserve the restriction of Q to Af^ if and only if 

trV'(x) = 0, 

for X € Mj-f. In particular, if the Hamiltonian h is a G-invariant function, then 
r2|jv/^ is an invariant volume form. 

3. Partial Lagrange-Routh Reductions 

Let (Q, I) be a natural mechanical system with Lagrangian I = ^{Kqq, q) — v{q), 
where the metric n is also regarded as a mapping n : TQ T*Q. The motion of 
the system is described by the Euler-Lagrange equations 

or by the Hamiltonian equations on the cotangent bundle T*Q with the Hamiltonian 
h{q,p) = i(p, K~^p) + v{q) being the Legendre transformation of I. 

Let G be a compact connected Lie group acting freely on Q and tt : Q — s- 
B — Q/G he the canonical projection. The G-action can be naturally extended to 
the Hamiltonian action on T*Q: g ■ {q,p) ~ {g ■ q, {dg^^)*p) with the momentum 
mapping $ given by (e.g, see (TH]) 

Here is the vector given by the action of one-parameter subgroup exp(t^) at q. 
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For Lagrangian systems, it is convenient to work with tangent bundle reductions. 
Let Vq = {^q I <^ G 0} be the tangent space to the fibber G ■ q [vertical space at q) 
and V ~ UqVq be the vertical distribution. Then 

(r*Q)o = $-1(0) = U,annV,, annV, = {p e T^Q \ (p,^,) = 0, ^ € 0}. 

Consider the horizontal distribution H — UqHq C TQ orthogonal to V with 
respect to the metric k. Equivalently, H is the zero level-set of the tangent bundle 
momentum mapping 

(25) n = <f^\0), ($,(g,(?),0= (1^,^,) =(^9<?,e?), eeS- 

Since Kq{Ti.q) = annVq, we see that Ti. is invariant with respect to the "twisted" 
G-action 

(26) goiq,X) = {g-q,K;.loidg-')*oKqiX)), XeTqQ, 

that is the pull-back of canonical symplectic G-action on T*Q via metric k: 

K 

TQ — > T*Q 

90 i I g- 

TQ — * T*Q 

K 

From Theorem [1] we obtain 

Theorem 3. (i) The horizontal distribution \25\) is an invariant submanifold of 
the Euler- Lagrange equations !124\ I if and only if the potential v and the restriction 
Ku of the metric k to Tl are G -invariant with respect to the action \26\) . 

(ii) // Ti. is an invariant submanifold of the system (Q, I) then the trajectories 
q(t) of the natural mechanical system with velocities q{t) that belong to Ti project 
to the trajectories b(t) = TT{q{t)) of the natural mechanical system (_B, L) with the 
potential V{TT{q)) = v(q) and the metric K obtained from k-^ via identification 
n/G^TB. 

Note that when k is G-invariant, the twisted G-action ([26]) coincides with usual 
G-action: g ■ (q,X) — {g ■ q,dg(X)) and the induced metric K is the submersion 
metric. In this case Theorem [3] is exactly the classical method of E. J. Routh for 
eliminating cyclic coordinates [27l [21 [21] . 

Definition 2. By the analogy with the Lagrange- Routh reduction we shall call 
the procedure of passing from the Lagrangian system (Q, I) to the system (B, L) a 
partial Lagrange- Routh reduction. 

Remark 5. The distribution (|25p . in general, is nonintegrable. One can interpret 
the restriction of the Lagrangian system (Q, I) to as a nonholonomic system 
(Q, I, Ti.) with a property that the reaction forces are equal to zero. 

3.1. Geodesic Flows on Compact Lie Groups. Let G be a compact Lie group, 
g be the Lie algebra of G and (•,•) be a Adc-invariant scalar product on g. Let 
a e g be an arbitrary element and let 6 belongs to the center of Qa- Here Qa = 
{ry e g, [a, 77] = 0} is the isotropy algebra of a. Let g = ga + be the orthogonal 
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decomposition. Consider the linear operator (so called sectional operators [13]) 
Aa,b,c : ^ 0, defined by 

Aa,6,c(0 = ad^^ oadfcoprJO + ^iP^sa 0> 

where prj, and pr^^ are the orthogonal (with respect to (•, •)) projections to and 
0a, respectively and C : 0a ^ fla is symmetric. We can always find b and C such 
that Aa.b,c is positive definite. 

Identify 0* with g by means of the scalar product (•, •) and consider the left- 
trivialization: 

r*G«zGx0 = {(5,e)}. 

Then the quadratic form 

(27) Vb,c = ^(AaAc(0,0 

can be regarded as the Hamiltonian of a left-invariant Ricmannian metric on G. 
Denote this metric by Ka,b,c- 

Let Ga be the adjoint isotropy group of the element a and consider the right 
Ga-action on G. With the above notation, the momentum mapping and its zero 
level-set are given by 

(28) $(g,0=pr,„(e), 

(29) (T*G)o«iGx5. 

Lemma 4. The metric K,a,b,c is invariant with respect to the Ga-action, i.e., the 
momentum map i2l^) is preserved along the geodesic flow 

(30) e= K,Aa,b,c(0], 9^g-KAc{0 
if and only if the quadratic form (^, C(^)) is Ada ^-invariant: 

(31) [e,C(O]=0, eesa. 

If a is regular element of the Lie algebra 0, then 0a is Abelian and the condition 
(|3T|) is satisfied for an arbitrary operator C. 

Suppose (|3ip holds. Then we can project the metric Ka,b,c to the homogeneous 
space G/Ga, that is to the adjoint orbit of a: G/Ga ~ G(a) = {Adg(a) | g E G}. 
Since we deal with the right action, the vertical distribution is left-invariant: Vg = 
g ■ Qa, while from the definition of Ha,b.Ci the horizontal distribution is Tig = g ■ d 
and the submersion metric docs not depend on C. 

Denote the submersion metric by Ka^b- The cotangent bundle T*0{a) can be 
realized as a submanifold of x 

(32) T*Oia) = {ix,p) I X = Ad<,(a),p £ 0^}, 

with the pairing between p G T*0{a) and rj g TxO{a) given by p{ri) = (p, 77). Let 
bx — Adg b. Then the Hamiltonian and the geodesic flow for the metric Ka,b in 
redundant variables {x,p) are given by (see [6]) 

(33) Ha,bix,p) = ^{a,db^p,a.dxp) = - ^{ad^ a.db^ p,p), 

(34) X ^ -ad^adb^p ^ [[bx,p],x], 

(35) p = -ad^^[p, [x, [bx,p]]] + pr [[bx,p],p]. 
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Now, let us perturb the metric Ka.b,c as follows. Take i5 = (B5(prg^),^) + 
hi^siWg^ 0>0> : c) ^ g^, Cs : Qa ^ 3a, such that 

(36) hsig,0^ha.b.c + S=^{AsC,0 



is positive definite. Then ([36]) will be the Hamiltonian function of the left-invariant 
metric that we shall denote by ks. Since hs — ft.|(T*G)o) the geodesic flow of ks 

(37) e=K,A,(0], <? = 5-A,(0. 

has the invariant relation (I29p . Hence we can perform the partial reduction. The 
metrics ks and Ka,b,c induce the same metric Ka,b on the orbit 0(a), but their 
horizontal distributions and H are different for B 7^ 0: 

In the case when a is a singular element of g we can take S = ^{Csipig^ ^), £,) such 
that d is not Adg^ -invariant. Then the perturbed metric has the same horizontal 
distribution as the non-perturbed one: Tig = g ■ I). 

The geodesic flow (p4|). ([35]) is completely integrable in the non-commutative 
sense. Moreover, the system is also integrable in the usual commutative sense 
by means of analytic functions, polynomial in momenta and which can be lifted 
to commuting Ga-invariant functions on T*G (see [5j [22l [6]). Thus, according 
Theorem [2] we obtain 

Corollary 5. The equations I130\) and fj*?] ) have the same invariant isotropic foli- 
ation of (2^. 



3.2. Local Description. To clear up the difference between geodesic flows ([30]) 
and ([57]) . let us write down the problem in local coordinates. 

Suppose a is a regular element of g. Then Ga ~ T" ia a maximal torus. Locally, 
we have T*G « T*0{a) x r*T". We take coordinates {q, ip,p, 0) in T*G such that 

{q,p) ^ {qi, ■ ■ ■ ,qni,Pi, ■ ■ ■ ,Pni) and ((p, 0) = ((pi,...,(/?„,0i,...,0„) 

are canonical local coordinates on T*0{a) and T*T'^ . Locally, the right T"-action 
is given by translations in coordinates and (r*G)o — {{q, fjP, 0)}. 
The Hamiltonians ha,b,c and kg are of the forms: 

where A*-', i?*", C"^ do not depend on the variables ipa- 
The partially reduced system 

* = ^=E^>^' P^ = -^-§^ = -T.^4^P^P^^ z = l,...,m 

is completely integrable and we can treat q{t) and p{t) as known functions of time. 
Here Ha.b — \ ^ A^^'PiPj is the reduced Hamiltonian ([33]) . 

The equations of the geodesic flows ([30]) and (|37p on (|29p in variables qi^pi are 
given by (|38p . The equations in variables (^a, respectively, are given by: 



= ^ B'"{q)pi and = " = 1' ■ 



. ,n. 
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While the first system is solvable, the second one, generically, is not. Globally, ([29]) 
is foliated on invariant tori with quasi-periodic and non-quasi-periodic flows of (|30|) 
and ([57)1 . respectively. 

4. n-DlMENSIONAL HeSS- ApPEL'ROT SYSTEM 

Consider the motion of an n-dimensional rigid body around a fixed point O = 
(0, 0, ... , 0) in the n-dimensional Euclidean vector space (R", (•, •)). The configura- 
tion space of the system is the Lie group SO{n): the element g e SO{n) maps the 
moving coordinate system (attached to the body) to the fixed one (e.g., see [HI). 
Let Fi, . . . , Fn, El, ... , En and /i, . . . , /„, ei, . . . , e„ be the orthonormal bases at- 
tached to the body and fixed in the space regarded in the space frame and moving 
frame: 

Ei^h = (l,0,...,0,0)^,...,S„ = /„ - (0,0,..., 0,1)^, 
El ^ g ■ ei, . . . ,En = g ■ e„, Fi = 5 • /i, . . . , F„ = 5 • /„. 

We can consider the components of the vectors ei, . . . , e„ (or Fi, . . . , Fn) as redun- 
dant coordinates on SO{n). 

For a path g{t) € SO{n), the angular velocity in the body frame and angular 
velocity in the space frame are defined by uj{t) = g^^ ■ g{t) e so(n) and Vl[t) = 
g ■ g~^ =^ g -uj ■ g~^ — Adg w, respectively. From the conditions ~ Ei ~ g ■ ei+ g ■ Ci 
and Fi — g ■ fi, the vectors ei, . . . , and Fi, . . . ,Fn satisfy Poisson equations 

ei = -Lo ■ a, Fi^n-Fi, 1 = 1,..., 71. 

The kinetic energy of a rigid body is a left-invariant quadratic form ^(luo,uj), 
where I : so{n) — > so{n) is a non-degenerate inertia operator and {X, Y) — 
— ^tr(Xy) denotes the Killing metric on so{n). For a "physical" rigid body, Itj 
has the form luj + loI, where J is a symmetric n x n matrix [12j . We will relax this 
condition, considering an arbitrary positive definite operator. 

Further, suppose that the body is placed in the homogeneous gravitational force 
field in the direction and the position of the center of mass of a rigid body is 
pfn. Then the potential is v = p&dJl {fm^n)i where © is the gravitational constant 
and is the mass of the body. Let m = I cj be the angular momentum in the body 
frame and A = I By the use of Killing metric we can identify so{n) and so(n)* . 
Then the Hamiltonian in the left-trivialization 

(39) T*SO{n) SO{n) x so{n) = {{g,m)} 
reads 

h = i(TO,ATO) +p637l(/„,e„) 

and the equations of the system, in redundant variables (ei, . . . , e„, m), take the 
form of the Euler-Poisson equations 

(40) rh = [m, ui] + p&M /„ A e„, w = Am 

(41) Ci ^ -uj ■ Ci, i = 1, . . . ,n. 



12 



BOZIDAR JOVANOVIC 



4.1. Invariant Relations. Consider the orthogonal, symmetric pair decomposi- 
tion J + O of the Lie algebra so{n): 

i ~ span{/iA/j, I < i < j < n — 1} = so(n— 1), d = span{/iA/,i, 1 < i < n— 1}. 
Then we can write h as 

h = i(m6,Atme) + (m{,Bmo) + ^(too, AoTOo) + p©9Jl (/„, e„), 

where Aj — pr^ oA o pr^, A^ — pr^ oA o pr^, B = pr^ oA o pr^, m{ = prj to, too = 
prg m. 

Let SO{n — 1) be the subgroup of SO{n) with the Lie algebra Consider the 
right SO{n — 1) action (rotations of a rigid body "around" the vector /„). The 
momentum mapping and its zero-level set are given by $ = TO{ — prj to and 

(42) {T*SO{n))o = $"^(0) {(g, to) | TO( = 0} = 50(n) x c) . 

From the relations [J,^] C 0, [0,0] C 6, the orthogonal projections of equations 
(|40|) to ? and c) are given by 

(43) TO{ = [TOt,A{TOt] -f [to(,Btoe,] -I- [m£,,Aor7io] -I- [toj, , B'^toc] , 

(44) Too = [toj, AgTOo] -I- [m{,B'^TO(.] + [too,A{TO(.] -f [toe,,Btoo] +p©OT/„ Ae„. 

Hence, (j42p is an invariant submanifold if and only if [TOj,,Aj|TOa] = 0, i.e., 
(too, AaTOj) is the Adso(„_i)-invariant quadratic form on i) . Since SO{n) / SO[n—l) 
is a rank one symmetric space, Adso(„_i)-invariant quadratic form on 5 is unique 
(up to multiplication by a constant). Thus, we come to the following proposition 

Lemma 6. The submanifold is an invariant set of the Euler-Poisson equations 
P^ ), if and only if 

(45) prj, oAo prg = aldo, 

where Id;, is the identity operator and a G M. 

Suppose (|i5|) holds. In coordinates to^ = (to, A fj), the invariant set is 
given by 

(46) niij =0, 1 < i < j < ?i - 1. 
For n = 3, after usual identification (M'^, x) = (so{n), [•,•]) 

(0 -TO3 TO2 \ 

TO3 —mi , 
— TOl y 

the operator A take the form ([6|) and the invariant relation (|46|) become ([7]). 

Therefore, it is natural to call the rigid body system (j40ll . (|4T|) . where (j45|) holds, 
a n- dimensional Hess-Appel'rot system. 

4.2. Right Trivialization and Reduction to 5""^. In order to describe the 
partially reduced system we consider the problem in the right-trivialization 

T*SOin) «^ SO{n) x so{n) = {(g,M)}, 

i.e., in the frame fixed in the space. Instead of (ei, . . . , e„) we use (Fi, . . . , F„) as 
redundant coordinates on SO{n). Then the sphere S"""^ = SO{n)/SO{n — 1) can 
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be identified with the positions of the vector F„: 
S0{n-1) — > SO{n) 

, 7r(Fi,...,F„) = i^„. 

The equations of motion (|40)) . (|4T|) in the space frame take the form 

(48) M = p©9JlF„AS„, 

Fi^Vl- Fi, I = 1, . . . ,rt, 

where the angular momentum M and velocity 17 in the space are related to the 
angular momentum m and velocity u in body coordinates via M — Adg m and 
D, = AdgOj, respectively. Therefore 

(49) r2 = Adg(w) = AdgoAo Adg-i(M). 

Since Fi A Fj — Adg{fi A /j), the invariant submanifold (f42|) in the right- 
trivialization reads 

(T*50(n))o ~r {{9,T^g), -Dg = Ad<,(c))}. 

Let JCg ~ Adg (t) and X = AT^: + Xx> be the decomposition of AT G so{n) with 
respect to so(n) = JCg + 'Dg. Note that the orthogonal projection to Vg given by 

(50) Xv^ = (X • Fn) A 

On (r*50(n))o we have w = Ao(m) + B(m) = a • m + B(to). Thus 

Q.^VLk.^+Q.v^, 17k, = AdgoBo Adg-i(M), f^p, = a • Af, 

Furthermore, from ([50| we have • F„ = ilp^ • F„. Therefore, on (r*50(n))o, 
the Hess-Appel'rot system form the closed system in variables (F„,fiug): 

(51) tlr^^^ap'&mFnSEn, 

Fn — flx>g ■ Fn, 

From the second equation we get 

(52) fl-D^ = F„ A F„, i.e., M = -F^ A F„ 

a 

and the first equation can be rewritten in the form: 

(53) K = -ap&m E„ + nFn . 

Here the Lagrange multiplier fi is determined from the condition (i^„, Fn) = 1. 
Therefore we obtain: 

Theorem 7. The motion of the center of the mass of the n- dimensional Hess- 
Appel'rot system |^5y m f/ie space frame is described by the spher- 
ical pendulum equations i5S\) . 

The same statement, for the classical Hess-Appel'rot system is proved by Zhu- 
kovski (see [^[5]). 

It is well known that the partially reduced system ([55]) is completely integrable 
in the non-commutative sense. The tangent bundle TS*"^^ is foliated by two- 
dimensional invariant manifolds that are level-sets of the energy and integrals 

{Fn A Fn, E, A Ej) = {Fn, Ei){Fn,Ej) - {Fn,Ej){Fn,E,), l<i<j<n~l. 
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Using ((52|) . we see that the lifting of these integrals to {T* SO{n))o are compo- 
nents of the momentum map of the left SO{n — l)-action: 

(54) Gij = a ■ (M, EihEj), \<i<j <n-l 

(note that (jM]) are integrals on the whole phase space T*SO{n) as well). Whence 

Corollary 8. The invariant set 114 of the Hess-Appel'rot system is almost every- 
where foliated by (2 -\- (n — l){n — 2)/ 2) -dimensional isotropic invariant manifolds, 
level sets of the Hamiltonian function and integrals |5^[ ). 

The above considerations motivate us for the following definition 

Definition 3. We shall say that a Hamiltonian system (fTT|) satisfies geometrical 
Hess-Appel 'rot conditions if it has an invariant relation Q and the partially reduced 
system P3|) is completely integrable. 

The geodesic flows on compact Lie groups presented in Section 3 provide exam- 
ples of systems that satisfy geometrical Hess-Appel'rot conditions. Note that the 
condition that the partially reduced system is integrable differs from the notion of 
restricted integrability given in [11] . 

4.3. Reduction to (so(n) x M")*. The system dM]), (gH) is always left SO{n-l)- 
invariant (rotations of a rigid body "around" the vector e„). Denote 7 = e„, 
r — pfn- The equations (|40|) together with the last Poisson equation, 

(55) m = [m, uj] -\- ©9Jt r A 7, 7 = — w -7, uj = Am, 

can be seen as a left SO{n — l)-reduction of ([^0]) . (|4ip . This is a Hamiltonian 
system on the dual space of Lie algebra e(n) = so{n) x R" with respect to the usual 
Lie- Poisson bracket. For n = 3, after identification (|47|) . the equations ([55]) get the 
familiar form ([T]). 

Now suppose the Hess-Appel'rot condition ((45|) is satisfied. Then the Euler- 
Poisson equations on the invariant set coordinately read: 



ri-l 

E 



UJ. 



ijmjri - p&Tl-fi, i = 1, . . . ,rt - 1, 



-1 



(56) 7,; = -a7„mi„ - ^ ^^in^' i = 1, . . . , n - 1, 



n-1 

in = 7jmj„, 
j=i 



where ujij = (Bm, fi A fj), 1 < i < j < n — 1. The equations (|56[) always have three 
integrals 




(57) -^3 = / (TOm7j - mjn-yiY, 



which correspond to integrals As we sow above, the variable 7„ = (e„,/„) = 
{En,Fn) satisfies the vertical component equation for the motion of the spherical 
pendulum and can be found by quadratures. 
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In general, equations ()56|) have no smooth invariant measure. The non-existence 
of a smooth invariant measure is the reflection of the non-solvability by quadratures 
of the classical system (IT|), (H)). Namely, by the Euler-Jacobi theorem (e.g., see 
[2], page 131), the invariant measure together with the integrals ([2]), (|4]) would 
implies solvability of the system by quadratures. 

Example 1. As an illustration, consider the case n — A and Hamiltonian 

h = 2('^i"^23 + azmjg + 03771^2) + 2 ("^14 + "*24 + "^34) + 
+bimi2m,ii + b2mi2m24 + 63'7ii2m34 + p69Jl7„. 

Then the equations in variables TO12, 7^23, "^13 are 

"112 = nT-i3'Ti23(a2 - fli) + mi2{bim24 - b2mii), 

"l23 = nil2'71l3(a3 - 0.2) + m.l2(^2TO34 - ^3"l24) + mi^ibimu + &2TO24 + 63TO34), 

"ii3 = "i23?ni2(ai - 03) + "ii2(6iTO34 - 53TO14) - m23(6imi4 -I- 62'7i24 + b^ims,^), 
while the equations in other variables, on the invariant set 

(58) mi2 = mi3 = m23 = 0, 
take the form 

mi4 = -m24(feimi4 -|- 52TO24 + &3ni34) - p(S9Jl7i, 
m24 = mi4(6imi4 + 62ni24 + &3"i34) - p©9Jl72, 
"^34 = -p©97t73, 

(59) 71 = -a74mi4 - 72(6imi4 -f- 62"i24 + fe3"i34), 

72 = -a74m24 + 7l(6imi4 + 62"l24 + fe3"l34), 

71 = -07477134, 

74 = a(7i?Tli4 + 7277724 + 73TO34)- 

Together with ([FT]) equations ([55)1 have a supplementary integral 

(60) T12 = mi4j2 - "12471: 

implying the foliation on 3-dimensional invariant manifolds of (j58p . 

In redundant coordinates (Fi, F2, ^3,^4, M) the integral on T*S'0(4) reads 

(61) T12 - (m,/iA/4)(e4,/2)-(r77,/2A/4)(e4,/i) 

= (Af, Fi A F4) (F4 , F2 ) - ( Af, F2 A F4) {Ei , Fi ) . 

It is clear that T12 is independent of integrals ((54|) . Whence the 5-dimensional 
invariant isotropic manifolds outlined in Corollary [5] are foliated on 4-dimensional 
invariant level-sets of (pTj) . 

Note that the system ([5^ has an invariant measure if and only if 5i = 62 = 0. 
Also, recall that for ai — a2 = a, 61 = 62 = 63 = 0, the Euler-Poisson 

equations represent the motion of a dynamically symmetric heavy rigid body, 4- 
dimensional version of the Lagrange top [3| . 

Therefore, if 5i = 62 = 63 = and a\ ^ a2 ^ 03 then the equations ([55]) 
coincides to the 4-dimensional Lagrange top equations restricted to the invariant 
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level-set ([58]) . although the systems are different on the whole phase space. They 
have additional integrals 

^13 = mi473 - m347i (= (M, Fi A F3){Ei, F3) - {M, F3 A Fi){E4, Fi)) , 

(62) J^23 = ^2473 - ^3472 (= (M, F^ A F4){E4, F3) - (M, F3 A ^4) (^^4, F2)) . 

Since the equations (j59p preserve the standard measure in variables toi4, m2i, 77134, 
71,72,73,74 and the invariant manifolds given by integrals ((57)) . ([60)) . ([62)1 are two- 
dimensional, by the Euler-Jacobi theorem, they are solvable by quadratures. On 
the other side, the Hamiltonian and functions ([54l) . (f6T|) . (f62|) provide 3-dimensional 
invariant foliation of 9-dimcnsional manifold (T*S'0(4))o. 

4.4. L-A Pair. The non-commutative integrability of the n-dimensional Lagrange 
top is proved by Beljaev [3] and the L-A pair is given by Reyman and Semenov-Tian- 
Shansky [26]. The L-A pair of the similar form can be used for the Hess-Appel'rot 
system on so{n) x M". Let m*, uj* , 7*, r* be the so(n + l)-matrixes given by 

where is the zero n x n matrix. 

Then, under the Hess-Appel'rot conditions (|45p and (|46l) . the Eulcr-Poisson 
equations ()55|) are equivalent to the matrix equation with a spectral parameter 
A 

(63) L{X) = [L{X),A{X)], 

where L(A) = j* + Am* + X'^^&Tlr* and A(A) = lo* + XeMr* . 

Indeed, the straightforward computations shows that the terms with A'' and A^ 
in ([63]) are equivalent to the equations (|55p . The left hand side term with A^ is 
identically equal to zero. It can be proved that the right hand side term with A^ 
vanish, on the invariant submanifold (j46|) . if and only if (|45|) holds. 

On the invariant set (f46|) the spectral curve is 

(64) p{X, fi) = det(L(A) - fild) = (-a*)""' (m"* + m'^(A) + QiXf) = 0, 

PiX)=T2 + -X'T, + X^(^y , Q{X)^XT3, 
a \ a / 

where Tx^T^^Tz are integrals ([57)1 . The curve ([M]) is reducible (for n > 3) and 
consists of the n — 3 copies of the rational curve /i = and the singular algebraic 
curve 

(65) F: + ii^P{X)^Q{Xf 

The spectral curves of the form (f65|) appear in the algebro-geometric analysis 
given in piJlfTT]. 

In particular the above consideration leads to the L-A pair of the classical system 
different from those obtained in [9]. 
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Theorem 9. Under the Hess-Appel'rot conditions (0), the Euler-Poisson 
equations {Ip are equivalent to the matrix equation W3\) . where 



I 



m = 



A{X) 



V 





-Am2 
-71 








Atoi 

-72 






Am2 
—Ami 


-73 - A^eOJl^ 


Xp^m 




73 



7i 

72 

x^<&m-^ 

0.2 





UJ2 



-A/9©2Jl 



Remark 6. Recall that the spectral curve for the classical system given in [9] is 
reducible and consists of the rational and the elliptic curve. In our case, due to the 
involution r : (A, /i) 1-^ (A, —p), the singular spectral curve (|65p is a two-covering of 
the regular (for the generic values of the integrals genus-3 hyperelliptic curve 

r' = r/r : + uP{X) + Q{Xf = 

(e.g., see Lemma 1 in [10]). In the afRne part, the point (0,0) € F is an ordinary 
double point as well as the ramification point of the covering. 

5. Hess-Appel'rot system on (so(4) x so(4))* 

Dragovic and Gajic defined a class of systems of the Hess-Appel'rot type by using 
the analytical and algebraic properties of the classical system [TT]. Subtle differ- 
ence between Dragovic-Gajic systems and those constructed here has to be studied 
carefully, but there is a big coherence between these two approaches. Remarkably, 
the so(4) X so(4) system given in [11] considered on the whole phase space T* 5*0(4) 
satisfies geometrical Hess-Appel'rot conditions as well. 

We start from the example of the geodesic flows given in Section 3. Suppose, in 
addition, the Ga-invariant potential force v{g) is given. Then remains to be an 
invariant manifold for natural mechanical systems with kinetic energies Ka,b,c a-nd 
K5. After Lagrange- Routh and partial Lagrange- Routh reductions, the systems 
project to the system with kinetic energy Ka.b and potential force V{x) on the 
adjoint orbit C(a), where V(x) = v(g), x = Adg a. 

The reduced system is a modification of ([M]) . (|35p by the potential force: 

(66) 
(67) 



x = - aAx adf,^ p = [[bx,p], : 
p^~a.d^^[p, [x,[ba:,p]]] 



dV{x) 
dx 



Here ^1 (x) , . . . , ^„ {x) is a base of Qx and Lagrange multipliers Xi are chosen such 
that trajectory {x{t),p{t)) belongs to (l32]) . For a proportional to h and a hnear 
potential V(x) = (cc, c), the system represents analog of spherical pendulum on 
adjoint orbit 0(a), which is integrable on an arbitrary orbit 0{a) (see Bolsinov and 
Jovanovic [7]). 

Now consider the above construction for the case G = 5*0(4), i.e, for the 4- 
dimensional rigid body motion about a fixed point (we use notion from the previous 
section). Let 



(68) 



6 = ( Ji + .h)a, a = 012/1 A /2 + 034/3 A f^. 
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We have so{4)a = span{/i A /2, /s A Ja} = so(2) ® so(2). Define 

C : so{A)a so{A)a, C(/i A/2) =2Ji/i A/2, C(/3 A /4) = 2J3/3 A /4. 
Then the left-invariant kinetic energy ha^b,c (see (P7|) ) takes the form 

J. /' N t2it2i'^1+'^3/2i 2, 2, 2\ 

Consider the perturbation of tlic metric 

(69) hs = ha.b.c + S = ^{Asm,m), 

S{g,m) = TOi2(Ji4"li4 - Jl3"723) + m-34(Jl3'77l4 - Jl4"^23) 

and the potential function 

V = (Adg-i (a), a) = (ai2ei A 62 + a34e3 A 64, a) 

(70) = (a,Adg(a)) = (a,ai2Fi AF2 +034^^3 AF4) 

(written in the coordinates (61,62,63,64) and {Fi, F2, F^, F4), respectively). 

The rigid body system with Hamiltonian h = hs + v, in the left-trivialization 
((39)) . takes the form 

(71) 771 = [777, A5(7T7)] + [ai26i A 62 + 03463 A 64, a], g = g ■ ks{m). 
The system (j71l) has two invariant relations 

(72) 77712 = 0, 77734 = 

and it is reducible to the oriented Grassmannian variety Gr+(4, 2) « 0{a) of ori- 
ented two-dimensional planes in the four-dimensional vector space. The diffeomor- 
phism ^E" : 0{a) Gr+(4, 2) is simply given by 

^'(a;) = Fi A F2, where x = Adg(a) = ai2J^i A F2 -I- 034^^3 A F4. 

Recall that ^1,^2,^3, F4 is the moving base regarded in the space frame, so Fi A-F2 
is the oriented two-plane attached to the body. 

From ([55]) . ([70)1 we get the reduced Hamiltonian 

Ji + J3 

H{x,p) = Ha,b{x,p) + V{x) = {[x,p], [x,p]) + {x,a) . 

After calculating the Lagrange multipliers in ((66| , (IGT)) the partially reduced system 
on r*Gr+(4,2) become 

(73) x^i.h + J3)[[x,p],xl 

(74) p== (Ji + J3)[[a;,p],p] -a + pr^^a. 

It follows from Theorem 4 that the reduced system is completely integrable by 
means of integrals that can be extended to commuting S'0(4)a-invariant functions 
on r*S'0(4). Hence we get the following qualitative behavior of the system 

Theorem 10. The rigid body system ^71\ l satisfies geometrical Hess-AppeVrot con- 
ditions: the partial reduction of the system from to the oriented Grassmannian 
variety Gr+(4, 2) is completely integrable pendulum type system jlS^ , J 7^1 ). Further, 
10- dimensional invariant manifold |y^[ ) is almost everywhere foliated by invariant 
6-dimensional Lagrangian tori that project to the ^-dimensional Liouville tori of 
the reduced system \73^ , |74[ ). We need two additional quadratures to solve the 
reconstruction problem. 
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By introducing 7 ~ a\ie\ A 62 + 03463 A 64, from ([7T|) we can write the closed 
system in variables (m,7): 

(75) m = [to,cj] + [7,0], 7=[7,w], = A5(m). 

Regarding 7 as the free so(4)-variable, the system (ffS]) becomes the Hamiltonian 
system with respect to the Lie-Poisson bracket on the dual space of the semi-direct 
product so(4) x so(4) (see Ratiu [25]). Moreover, with the above choice of we 
have 

/ Ji Ji3 \ 



(to) = Jto + m J, J 



Ji J24 
^13 J3 

V J24 J3 y 

and ((75|) coincides with the Hess-Appel'rot system defined in [TT] , 



Remark 7. A similar statement can be proved for the Hess-Appel'rot systems on 
so(n) X so(ji) for n > 4 (see |TTI), where one should take 

6 = (Ji + J3)a, a = 012/1 A /2. 

Then so(n)a = so(2) so(n — 2) and the systems are partially reducible to the 
orbits 0(a), which are now diffeomorphic to the oriented Grassmannian varieties 
Gr+(n, 2) of oriented two-dimensional planes in the n-dimensional vector space. 



Relations (|72p define 6-dimensional invariant manifolds within generic symplectic 
leafs of (so(4) x so(4))*. Appart of the Hamiltonian function, the system ([75]) has 
a supplementary integral (see [Tl) ) 

T = 734012 + 712034 + {J\ + J3)('7ll2m34 -|- m23mi4 - 7711377124) 

implying the foliation on 4-dimensional invariant varieties. The L-A pair represen- 
tation as well as the algebro-geometric integration of ([75]) up to two quadratures 
are given in [llj . 
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